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Abstract 

We investigate the dynamical properties of a system of interacting magnetic dipoles disposed 
in sites of an sc lattice and forming a cubic-shaped sample of size determined by the cube edge 
length (N — l)a (a being the lattice constant, N representing the number of dipolar planes). The 
dipolar field resulting from the dipole-dipole interactions is calculated numerically in points of the 
axis connecting opposite cube face centers (central axis) by collecting individual contributions to 
this field coming from each of the N atomic planes perpendicular to the central axis. The applied 
magnetic field is assumed to be oriented along the central axis, magnetizing uniformly the whole 
sample, all the dipoles being aligned parallelly in the direction of the applied field. The frequency 
spectrum of magnetostatic waves propagating in the direction of the applied field is found numer- 
ically by solving the Landau-Lifshitz equation of motion including the local (nonhomogeneous) 
dipolar field component; the mode amplitude spatial distributions (mode profiles) are depicted as 
well. It is found that only the two energetically highest modes have bulk- extended character. All 
the remaining modes are of localized nature; more precisely, the modes forming the lower part of 
the spectrum are localized in the subsurface region, while the upper-spectrum modes are localized 
around the sample center. We show that the mode localization regions narrow down as the cube 
size, N, increases (we investigated the range of N=21 to iV=101), and in sufficiently large cubes 
one obtains practically only center-localized and surface-localized magnetostatic modes. 

PACS numbers: 75.30.Ds, 75.40.Gb, 75.75.+a 

Keywords: magnetostatic modes; magnetic nanograms; mode localization; demagnetizing field 
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I. INTRODUCTION 



The development of nanotechnology in recent years allowed to design nanometric ferro- 
magnetic materials of any required shape with a large size variation ranged. Small magnets 
and particles have raised an increasing interest due to their potential application in magnetic 
random access memory (MRAM) elements and ultrahigh-density storage. This interest also 
involves the sample ability to rapidly change the magnetization state in order to optimize the 
individual recording time. Particularly, magnetostatic modes are of practical importance in 
this matter (especially in magnetization reversal processes), by reason of their low frequency 
range. 

The small magnetic elements used in signal processing devices are of nanoscale size 
and mostly non-ellipsoidal shape. In such small magnetic particles the internal field non- 
uniformity due to the particle shape has to be taken into account. Of course, the dominant 
mechanism in the creation of such a nonhomogeneity are long-range dipolar interactions. 
The dipolar size and shape effects present in the magnetostatic wave spectra of small nonel- 
lipsoidal particles result in additional power losses^ very important in device applications. 

The nonuniform demagnetizing field in a normally magnetized disk was studied long ago 
by Yukawa and Abe^, with the assumption of uniform saturation magnetization throughout 
the disk. Analyzing the ferromagnetic resonance (FMR) spectrum of the magnetostatic 
waves in a normally magnetized YIG disk, the authors found that these waves were excited 
at the disk edge and propagated radially towards the disk center with increasing wave 
numbers. 

The dynamical properties of small magnetic elements have been actively investigated 
recently in studies using different methods and approaches. Resonance methods combined 
with surface tunneling microscopy techniques enabled several team o 4 i 5 i 6 to observe magne- 
tostatic and exchange modes on a very local scale; time-resolved and space-resolved MOKE 
experiments^ were used to study the local magnetization dynamics as well. These differ- 
ent methods give good indications of the localized magnetostatic mode existence, especially 
when the applied field is strong enough to induce a magnetization saturation throughout 
the sample. 

The localized magnetostatic mode dynamics has been investigated by several authors. 
Berkov et a/. 8 studied the spin wave frequency spectrum and spatial distribution in thin 
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/xm-sized magnetic film samples. The lowest-frequency modes were found to correspond to 
oscillations restricted to the boundary regions, which start to 'propagate' inside a rectangular 
magnetic sample as the frequency increases. Tamaru et a/.— applied a spatially resolved FMR 
technique to quantized magnetostatic mode imaging in small magnetic structures. They 
identified quantized magnetostatic modes in the obtained spectra, and from the observed 
spatial distribution of magnetization response at each mode peak deduced that the number 
of mode nodes decreased with increasing bias field. However, in their data interpretation 
some discrepancies were found between the measured mode frequency values and those 
calculated on the basis of the Damon-Eshbach theory^. An exact theory of spin wave mode 
quantization in small magnetic structures needs to be developed. 

The dipole-exchange spectrum of the discrete spin wave modes resulting from the bias 
magnetic field inhomogeneity inside a nonellipsoidal rectangular magnetic sample has been 
studied recently by Guslienko et al. n . The authors found that the strong inhomogeneity of 
the internal bias magnetic field along the magnetization direction results in spin wave mode 
localization at either the edges (exchange localization) or the center (dipolar localization) of 
the sample. We will show in this paper that in a cubic-shaped sample with an internal field 
nonhomogeneity resulting solely from pure dipolar interactions these two types of magne- 
tostatic mode localization are also present (it means magnetostatic modes are found to be 
localized not only at the cube center, but also on the cube surfaces). 

Our theoretical approach is based on a modified Draaisma-de Jonge treatment used in 
the dipolar energy calculation. In the original Draaisma-de Jonge treatment a ferromagnetic 
film is considered as a set of discrete magnetic dipoles regularly arranged in a crystalline 
lattice. The dipolar energy is calculated by collecting contributions from each dipolar lattice 
plane assumed to be parallel to the film surface. The dipoles within each plane are divided 
into two sets: those within a circle of radius R and those beyond that circle; the total 
contribution of the former set is calculated discretely through summing over the dipoles, 
while the contribution of the latter set is evaluated by integration (the circle radius should 
be large enough to assure reliable final results). In the approach used in our study, the 
discrete summation is performed over all the dipoles within a given dipolar plane; thus, the 
integration is avoided, and no approximation is involved in our dipolar energy evaluation. 

The reason to use a discrete approach in the dipolar field calculation is our general 
belief^ that the continuous models^ used in deriving magnetostatic modes become rather 



4 



questionable when applied to small samples containing a finite number of elementary spins. 
The complexity of introducing the long-range spin-spin interactions into the discrete ap- 
proach is balanced by the restriction of our study only to the case of a uniform mode within 
each dipolar plane, which reduces the 3D problem to a single dimension (this assumption is 
experimentally justified by the observed weak in-plane variability of the magnetic modes^). 



II. A FINITE SYSTEM OF MAGNETIC DIPOLES IN PLANAR ARRANGE- 
MENT 

We shall consider a system of magnetic moments n? arranged regularly in sites r of a 
simple cubic crystal lattice. The system is assumed to form a rectangular prism with square 
base (see Fig. [Qt). Let the prism base determine the (x, y) -plane of a Cartesian reference 
system, with the z-axis perpendicular to this plane. The reference point (0, 0, 0) shall be 
placed in the central site of the prism bottom. 

Let us calculate magnetic field kg "produced" by all the prism dipoles in a site indi- 
cated by internal vector R. According to the classical formula (obtained using the linear 
approximation^), field can be expressed as follows (in the SI units): 

^12 



Z 1 3(r-fl)(/2Hr-fl))-/ZHr-fl| 
* 4vr |f_ R\5 ' ^ ' 
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the above sum involving all the sites except the reference point (i.e. the site with position 
vector f=R). The lattice planes parallel to the prism base shall be numbered with index 
n G (0, N — 1) (see Fig. Hh), and the sites within each plane indexed with vector r\\ = 
a[pi + qj], defined by integers p,q G (—L,L). This means that the position of the sites in 
which the magnetic moments are located, indicated by vector r, shall be defined by a set of 
three integers, (p,q,n): 

r = [r\\,an] = a\p,q,n], p,qe(—L,+L) and n G (0, N — 1), (2) 

a denoting the lattice constant. Thus, the considered prism contains N(2L + l) 2 magnetic 
moments. Below we shall focus on the magnetic field on the z-axis only, assuming its 
direction to be solely allowed for magnetic wave propagation. Hence, we put R = a[0, 0, n'], 
where n' G (0, N — 1), and re-index the dipole field: h n i = hg. By vector decomposition: 

f — R = a [pi + qj + (n — n')k^ , (3) 
5 



equation JTJ becomes : 
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3a(pi + qj + (n — n')k) (jl? ■ a(pi + qj + (n — n')kfj 



{ a \/P 2 + q 2 + (n — n') 2 ) 5 

jl?a 2 (p 2 + q 2 + (n — n') 2 ) 



aJp 2 + q 2 + (n — n 



An 



E' 



+ /i|pg + /^p(n - n')) + j(p?pq + /ifq + ^pq(n - n')) 
a 3 {p 2 + q 2 + (n — n') 2 )l 
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(4) 



the triple sum over indices p, q, n I X) = J2 replacing the single sum over vectors r; the 
summing range is as defined in ©, and the prime (') at the sum symbol means that the ref- 
erence point [0, 0, n'] is excluded from the sum. As the site distribution within each plane n 
is symmetric with respect to site (0, 0, n), each site (p, q) has its counterpart (— p, —q). Con- 
sequently, in sum (pEJ), all the terms in which p and q appear in odd powers are compensated, 
and (H|) becomes: 

^ip%p 2 + jfJ^q 2 + k[Lf{n — n') 2 — ip%p 2 + q 2 + (n — n') 2 ) 



p,q,n 



a 3 (p 2 + q 2 + (n — n') 2 )^ 
jf4(p 2 + q 2 + (n- n') 2 ) + kfip(p 2 + q 2 + (n - n') 2 ) 
a 3 (p 2 + q 2 + (n — n') 2 )2 
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ifi'p(2p 2 — q 2 — (n — n') 2 ) + jfMp(—p + 2g 2 — (n — n') 2 ) 
a 3 (p 2 + g 2 + {n — n') 2 )% 

^ kfxp(-p 2 -q 2 + 2(n - n'f ) 



(5) 

a 3 (p 2 + q 2 + (n — n') 2 )^ 
Further calculations shall be performed using double sums over p and q, defined as follows: 

(n — n') 2 
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p 2 + q 2 
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(7) 



the equality of the two first sums in (J7J) is a consequence of the fact that the z-axis is a 
fourfold symmetry axis for each lattice plane (see Fig. Additionally, we shall assume 
that all the magnetic moments within a single plane n are identical, i.e.: 

fin = fi[ P ,q,n], for any p and q. (8) 
By including relations (JSHEj) into (0), magnetic field h n > can be expressed as follows: 



T7 



(9) 



Note that by assuming /2 n = fi? all the magnetic excitations propagating in plane (x, y) are 
excluded from our analysis, and z axis becomes the only direction of propagation allowed. 
In order to obtain a simpler expression of h n i we introduce a symmetric matrix whose 



elements, D nn i, are defined as follows: 



= _ _ l(p 2 + q 2 )-(n-n') 2 

J-'n,n' — <Jn,n' *n,n' / , k 



p,g \j)2 j rq 2j r (j l _ n iy 



(10) 



With this matrix, the magnetic field reads: 



D 



W n + hi - 2AK 



n,n' 



(11) 



It is convenient to introduce here the notion of magnetization, a phenomenological quantity, 
which in the considered case of simple cubic lattice can be defined as follows: 

M n = fi n / 'a 3 . 



Then, (JTTJ) becomes: 



47T 



(12) 
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It should be remembered that site (0, 0, n') is excluded from the sums appearing in the 
equation ([TO]): to avoid any ambiguity, this shall be noted explicitly by rewriting (|T3|) in its 
expanded form: 
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iM*+jM«-2kM*, 



(14) 



The second term on the right of Eq. (JUjl has been calculated for n = n', and thus site 
fu = (0,0) should be excluded from the sum over p and q (as indicated by the prime at the 
sum symbol). Equation (jTUjl . defining elements D n>n /, can be interpreted as the definition 
of a dipolar matrix D composed of these elements; as we shall see this matrix shall play an 
important role in deducing properties of magnetic modes. As a consequence of (|1U|) matrix 
D is symmetric, i.e. D n , n i = D n ^ n . Let us introduce a new variable, defined as follows: 

5 = n-n'; (15) 

5 measures the distance between planes n and n' . Consequently, we can write: 

D n>n , = D 5 = D n , in = D„ s ; (16) 

1 (V , 2\ _ r2 

m [p 2 + g 2 + 5 2 ] 5 

Note that the considered system consists of a finite number of planes (index n taking 
values n = 0, 1, 2, • • ■ , iV — 1), so the set of values available to 5 depends on the reference 
plane, n', with respect to which the distance is measured. However, the following condition 
must always be satisfied: 

<n + 5 < N -1. (18) 
A sum over all the system planes shall be in use below; this sum shall be denoted as: 

E n ', (19) 

s 

superscript n' indicating that the summing is performed on the planes neighbouring with 
n', including plane n' . Therefore, S takes the following values: 

5 = 0, ±1, ±2 ... (20) 

its lower and upper limits being determined by condition (fT8*j) . 

Up to now, the direction of dipole arrangement has not had much importance in our 
reasoning. Now we shall consider the case with dipoles arranged along the z-axis only. 

III. MAGNETIC DIPOLES ALIGNED ALONG THE DIRECTION OF MAGNE- 
TO STATIC MODE PROPAGATION 

In this paragraph we shall consider a magnetic prism placed in a static magnetic field, 
H Q , applied along the z-axis (Fig. [TJ. Field H is assumed to be strong enough to arrange 
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all the magnetic moments along the z-axis. Then, the magnetization vector can be regarded 
as a superposition of two components (Fig. [TJ): static (parallel to the z-axis) and dynamic 
(lying in the (x, y)-plane): 

M n = M s k + rh n ; (21) 

M s is the static magnetization, assumed to be homogeneous throughout the sample, and 
vector fh denotes the dynamic magnetization, perpendicular to Ms- Similarly, the dipole 
field, h n ', can be resolved into two components: static, h s n , (parallel to the z-axis), and 
dynamic, h% (lying in the (x, y)-plane): 

h n , = h s n , + h d n ,. (22) 

These two components of the dipole field can be easily found from (JT3J). By replacing the 
third component of the magnetization vector with the static magnetization (i.e. by putting 
M* = Ms), and the two other components, M^ and M%, with the respective components of 
the dynamic magnetization, and m y n) the following formulae are obtained: 



n n , 



M s k, (23) 



hi' = ^E^n'm„, (24) 

element D n ^ n i being defined by fTUj) . 

The magnetic moment dynamics is described by the phenomenological Landau-Lifshitz 
equation (LL): 

^ = 1^M A x H eff Ji , (25) 

eff R denoting the effective magnetic field acting on the magnetic moment in site R. This 
effective field is a superposition of two terms only: the applied field, H , and the field h n i 
produced by the magnetic dipole system: 

H e ff,R = Heff,n' — H + h n i. (26) 

Considering (|22j) . we can write further: 

H effin , = (H + h s n ,)k + h d nl ; (27) 
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the above-introduced dipole field components (static and dynamic) being defined by 
and (Eg). 

The LL equation becomes: 



dm n t 



7/4) 



(M s k + 



m„' x 



H + h s n ,) k + h d n ) . 



(28) 



We shall solve it using the linear approximation, i.e. neglecting all the terms with m 
squared. Assuming the standard harmonic time-dependence of the solutions: fh n i ~ e~ tuJt , 
becomes: 



ium n i = 7/i 



M s k x h% + rn n > x (H + h s n ,)k 



or, using the properties of vector product: 



ium n i = 'JHok x Msh A n , — m n i (H + h 8 n ,) 



(29) 



(30) 



Through replacing the dipole field dynamic and static components with their explicit ex- 
pressions (J23|l and (jUJ) we obtain: 



iujm ri i = 7/ioA; x 



M s -y~ Y\ D n ym n - m n > 

47T - 



or, after bilateral multiplication by — 47r(7/i Ms) 1 : 



iQm„> = k x 



fin -2 -]Tl> 

n,rr 

\ n / n 



(32) 



and denoting the reduced frequency and the reduced field, respectively, defined as 
follows: 



n = 



Attlu 



and Qjj 



4irH n 



IVoMs M s 
Two complex variables are now introduced for convenience: 



m n = m n ± im l\ 



(33) 



(34) 



with these new variables, ()32|) splits into two independent identical scalar equations for m + 
and m~ ; this means we are dealing with magnetostatic waves polarized circularly. Therefore, 
it is enough to consider only one of these two equations, e.g. that for m + : 



Q m + = 777+ lQ H - 2 D n,n> ) ~ D n,n>m A % 
\ n / n 



(35) 
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The above equation can be rewritten as follows: 

fim+, = 777+, j Q H - 2 E D n,n> ~ D n ',n' I ~ E D n,n>m+ . (36) 

Matrix elements D n>n r mean elements of the dipolar matrix discussed in Section [HI with 
notations introduced there, the equation (JSfiJl becomes: 

Q m + = m+,Q n/ - E D & m l>±5i (37) 

(5=1,2... 

where we introduced the following abbreviation denoting the local field: 

Q n , = n H -D -2j2 n 'D6- (38) 

5=0 

It is also convenient to introduce at this stage the reduced static demagnetizing field - on the 
analogy of the reduced external field (f3l3|) - by using the definition (J2BJ): 

4-7T / 

= TrK> = -2 E Am* = -2 E " A- (39) 

M ^ n 5=0 

This allows us to rewrite Eq. (J38|) in the form: 

n n , = n H + n d + n s n ,, (40) 

with n d = -D and SI', meaning the (reduced) contributions to the local field coming, 
respectively, from the dynamical and statical parts of the demagnetizing field. Note that, 
according to the assumption made at the beginning of this paragraph, the eigenvalues Q 
(being reduced frequencies) correspond to magnetostatic waves propagating in the direction 
of the applied field, i.e. along the central axis shown in Fig. |21 

In the remaining part of our work we will be considering only the cubic-shaped samples, 
i.e. starting from this point we always assume 2L = N — 1. We also assume particular 
values for /io-f^o — 0.2T and Ms = 0.139 • lO^Am' 1 (YIG magnetization) with resulting 
value for the reduced field Qh = 14.374. However, we have to emphasis that selection of 
this particular value for fin is not essential for results to be presented in subsequent sections 
of this work, since the distribution of eigenvalues Q and profiles of modes associated with 
them are not sensitive to the choice of particular Qh value: the particular value of Qh only 
sets the whole spectrum in a given frequency region and if Qh changes the whole spectrum 
is shifted to another region, but the relative distribution of mode eigenfrequencies remains 
unchanged. 
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IV. MODE FREQUENCIES AND AMPLITUDE PROFILES ALONG THE CUBE 
CENTRAL AXIS 



We shall investigate magnetostatic excitations in a cube of size 40a (a being the lattice 
constant). Fig. |2] shows this cube in a reference system rotated by 90° with respect to 
that indicated in Fig. ^ (the z-axis, along which the atomic planes are counted, is now 
horizontal). The cube consists of 41 planes normal to the z-axis and numbered with index 
n, ranging from n=0 (the left face) to n=40 (the right face). The effective dipole field 
calculation procedure applied in the previous paragraphs allows to find the field in z-axis 
points only, i.e. along the cube central axis, passing through opposite cube face centers; this 
is the idea of the approximation used throughout this study, and henceforth referred to as 
central-axis approximation (CAA). 

The whole sample is assumed to be magnetized uniformly (the corresponding magneti- 
zation value being Ms) and in the direction of the external field, applied along the z-axis. 
The z-axis direction shall be also the only one allowed for propagation of the magnetostatic 
waves studied in this paper. With these assumptions, the problem of motion - to be solved 
on the basis of (|25|) - reduces to a single dimension in the space of variable n; the domain 
of the investigated motion is the interval n G (0, N — 1), between two opposite cube face 
centers (Fig. |21 shows this interval in close-up). 

Fig. El presents the discrete spectra of numerically calculated magnetostatic mode fre- 
quencies in a cube of variable size; the spectrum evolution with increasing N, or cube size, 
is visualized by the depicted frequency branches, each corresponding to one mode of a fixed 
number m. The plot shows clearly that only in the lowest N value range (N <50) the 
frequency spectrum changes in a significant way; above this range, the frequency values 
stabilize at levels independent of N. A striking feature is that the frequencies of the two 
highest modes, m = N—l and m = N, as well as those of the two lowest ones, m = 1 and m 
= 2, are pronouncedly separated from the rather uniform ,,band" formed by the other mode 
frequencies. In Fig. EJ showing all mode profiles in a 40a x 40a x 40a cube, these "detached" 
modes reveal quite distinct amplitude distributions, differing from those of the other modes: 
modes m = N - 1 and m = N appear to be of the bulk- extended {BE) type (with anti- 
symmetrical and symmetrical amplitude distribution, respectively), whereas modes m = 1 
and m = 2 are surface-localized (SL), the lower one being antisymmetrical, and the higher 



12 



one symmetrical. Besides, note that the amplitude arrangement is ferromagnetic in the BE 
modes, and antiferromagnetic in the SL modes. All the other modes, within the ,,band", can 
be qualified as localized (L), their maximum amplitudes localizing in some specific regions 
inside the sample; as these localization regions are found to vary with the mode number, 
the mode localization appears to depend on the mode frequency. 

To examine this relation in detail, let's note that the band modes can be divided into 
two groups: modes m = 3 4- 24 and m = 25 -4- 39, showing different localization regions. 
A characteristic feature of the first-group modes is a zone of zeroing amplitudes around 
the sample center; two non-zero amplitude regions are present at both sides of the central 
"dead" zone. Each of these non-zero regions can be divided further into two sub-regions: an 
outer one, situated at the border, with ferromagnetic arrangement of the amplitudes which 
fade towards the sample borders, and an inner sub-region, with mostly antiferromagnetic 
arrangement of the amplitudes, reaching a maximum at a certain point. It is apparent 
from Fig. HJd that the "ferromagnetic tail" extends as the mode energy increases, shifting 
the neighbouring localization region inside the sample; finally the localization reaches the 
sample center, the modes entering the other group (i.e. m = 25 4- 39). Typical for this 
group, the central localization region tightens around the sample center as mode energy 
increases. With respect to the above-discussed localization properties in both groups, the 
second-group modes can be qualified as center-localized (CL), and the first-group ones as 
non-central bulk-localized (NCBL) (or, alternatively, empty-center bulk localized). 

V. TWO TYPES OF MODE LOCALIZATION: SURFACE-LOCALIZED AND 
BULK-LOCALIZED MODES 

Let's study in detail the NCBL mode properties. An extremely interesting feature of this 
group is both amplitude and energy double degeneration, appearing simultaneously. This 
effect is illustrated in Fig. showing the mode profiles divided into two groups: symmetrical 
modes (S) are separated from antisymmetrical ones (AS) (for simplicity reasons, only \m + \ 
profiles are depicted). Note that the profiles of consecutive mode pairs: m—1 (AS) and 
m=2 (S), m=3 (AS) and m=4 (S), and so on are identical. This amplitude degeneration in 
AS-S mode pairs takes place throughout the whole NCBL group, to pair m=23;24 inclusive. 
Its disappearance above this pair is easy to explain: as by definition, in antisymmetrical 
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modes amplitude vanishes in the sample center, a necessary condition for their symmetrical 
counterparts to have identical amplitude distribution is a similar amplitude zeroing in the 
center; but this occurs in NCBL modes only. The amplitude degeneration is associated with 
that of energy levels; a double energy level degeneration in the NCBL group is apparent in 
Fig. |U1 (When studying Fig. |Ul note once more that the two highest modes, m = N-l and 
especially m = N, are apparently separated from the other mode energy levels.) 

Getting back to Fig. note that all the modes have two ferromagnetic tails et their 
borders; the extent of these tails depends on mode frequency. In both the symmetrical 
and the antisymmetrical mode groups, the highest-energy mode consists of two complete 
ferromagnetic tails, each extending from a sample border to the sample center. In the 
highest-energy symmetrical mode, m — N, the tails meet in the center with their maximum 
amplitudes; in the counterpart antisymmetrical mode, m = N-l, the tails meet in the 
center with their zero amplitudes. An interesting thing to note is that each of these two 
highest modes can be regarded as the ,, source" of a mode family formed by all the other 
modes within the given symmetry group. The generation of the lower modes from the 
initial BE modes, m = N or m = N-l, is based on ferromagnetic tail shrinking: receding 
from the sample center, the tails give way to the formation of a new mode segment, with 
amplitudes localized in the center (CL modes). The central localization zone broadens as 
energy decreases, to split into two separate (right and left) localization sub-zones, which 
migrate towards the borders following the continuously shrinking ferromagnetic tails. No 
more non-zero amplitudes appear in the vacated center, which therefore remains "dead"; 
this corresponds to the formation of NCBL modes. When the non-zero amplitudes reach 
the borders (which occurs in the lowest-frequency modes), the entire inside is dead, and the 
corresponding modes are of surface-localized nature. 



VI. SIZE-INDUCED LOCALIZATION EFFECT 



Let's increase the cube size to examine its effect on the mode energy spectrum and 
profiles. Fig. IHb shows a discrete mode frequency distribution calculated for a cube of size 
2L=100; this distribution can be compared to that obtained for a cube of size 2L=40, shown 
in Fig. EK- Besides the obvious increase in the density of states (resulting from the sample 
enlargement), the larger cube retains the basic characteristics observed in the smaller cube, 
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namely: the double degeneration in the lower-spectrum modes, its absence in the upper 
modes, and the outstanding of the two highest modes. Note that also the mode frequency 
range remains virtually unchanged, and that with the increased density of states, the linear 
character of the mode frequency vs. mode number dependence is more apparent. 

However, this "invariance" of the energy spectrum does not translate into the amplitude 
profiles. Fig. [7| presents a juxtaposition of profiles of the five lowest modes and the six 
highest ones in cubes of two different sizes: 2L=40 and 2L=100. The profile invariance 
is found to be limited to the BE modes only (m=l;2) whose profiles remain unchanged in 
spite of the size increase. The other states reveal apparently changed localization: the NCBL 
modes in the 2L=100 cube, showing substantially larger dead regions (with respect to the 
2L=40 cube), are practically of sub-surface localized (SSL) nature; in the CL modes, the 
localization zone has tightened around the strict sample center. Thus, we can anticipate that 
further size increase shall result in a deepening of the above-mentioned localization changes, 
and that consequently, the magnetostatic mode spectrum in large cubic samples shall consist 
of three clearly different mode groups: (a) two highest-frequency bulk- extended modes (AS 
and S), (b) center-localized modes with intermediate frequencies, and (c) lowest-frequency 
surface-localized modes. Another way of illustration of this size-induced localization effect is 
given on the Fig. |HJ 

VII. INTERNAL DIPOLAR FIELD NONHOMOGENEITY AS THE SOURCE OF 
MODE LOCALIZATION 

Obviously, we would like to understand why almost all of the modes forming the magne- 
tostatic spectrum (except for the two highest ones) are of localized nature. We can suppose 
that this property is due to the long-range character of dipole interactions, as the latter as 
known to produce spatially nonhomogeneous demagnetization field in bounded non-ellipsoid 
systems^. In fact, this is the case of the cube considered in this paper: plotted with bold line 
in Fig. El the local field Q n , given by Eq.(|3S|). shows nonhomogeneous spatial distribution 
(along the central axis), reaching a maximum in the sample center, and fading symmetrically 
towards the borders. The magnetostatic mode profiles, plotted in the same energy scale, 
have been superimposed on this local field plot. Strikingly, the shape of the suppressed 
amplitude area in the bulk resembles exactly that of the nonhomogeneous distribution of 
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Q n , and all the NCBL mode frequencies lie within the local field value range. The CL mode 
frequencies coincide with the top of the Q n curve, i.e. the region where the local field gra- 
dient is highest. The BE mode frequencies lie above the Q n value range, and one can think 
that this is why they are hardly sensitive to the local field variations. 

VIII. FINAL REMARKS AND OUTLOOKS 

In their recently performed systematic studies of FMR in patterned submicron rectangular 
permalloy films, Zhai et al^ observed multipeak spectra interpreted as spin-wave resonance 
due to the non-uniform magnetization of the studied samples. We expect a similar study 
can be performed on cubic structures as well, with the aim of revealing different types of 
the above- described localized modes. 

One of the practical aspects of our study concerns the size- independence of magnetostatic 
mode frequencies found in cubic-shaped samples. This is in contrast with the recent report^, 
where FMR spectra were measured for a series of single crystal yttrium-iron-garnet (YIG) 
films of thickness ranging from 35 /im to 220 /*m. The spectra were interpreted numerically 
on the basis of a model similar to ours, i.e. taking into account the strongly inhomoge- 
neous internal field, but the resonance fields (measured for both in-plane and perpendicular 
configurations) were found to be sensitive to the film thickness. This effect is, in fact, not 
predicted by the results of our study, performed on cubic-shaped samples: almost identical 
local field fl n profiles (depicted with bold lines in Fig. EJ) were found inside two cubes of dif- 
ferent sizes, 2L=40 and 2L=100, and moreover, the mode frequency spectrum distribution 
with respect to Q n is almost identical in both cases. Therefore, we may conclude that the 
above-discussed size-independence of our magnetostatic mode spectrum resulted from the 
assumed cubic shape of the studied sample. 

A problem to be solved next is the mode localization behaviour when the sample is 
,, deformed", losing its cubic symmetry. In our next paper we are going to calculate the 
magnetostatic mode spectrum in a rectangular finite thin platelet. Our preliminary results 
indicate that the magnetostatic mode localization effects described here for cubes are present 
in rectangular platelets as well. We believe, in fact, that this effect is just underlying the 
so-called 'bimodal' statistical distribution of the internal field, recently revealed in finite 
YIG thin films by Pardavi-Horvath and Yan 2 , and found to be clustered around two major 
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components - a low-field peak corresponding to the central (volume) part of the sample, 
and a higher internal field related to the surface - which would correspond exactly to CL 
and SL modes disscussed in this work. Closer elucidation of this analogy will be the subject 
of our subsequent paper. 

Another question requiring future investigation is the mode localization dependence on 
mode propagation direction. In another study we are going to examine magnetostatic modes 
propagating perpendicularly to the applied field. 
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FIG. 1: (a) The prism sample considered here; the case in which the applied field, Hq, sets the 
magnetic moments in the direction perpendicular to the prism base (i.e. along the z-axis). The 
prism 'thickness' is (N—l)a, and the square base side width is 2La (a denoting the lattice constant). 
The magnetostatic waves are assumed to propagate along the z-axis (i.e. in the direction of the 
applied field), (b) The planar prism model used in our calculations. 

FIG. 2: We consider a cube of edge size 40a, magnetized uniformly along the z-axis by the applied 
magnetic field (Ms being the magnetization). Magnetostatic modes are assumed to propagate only 
in the direction of the applied field. The dipolar field entering the equation of motion is calculated 
numerically along the central axis (brought into close-up in the graph). 

FIG. 4: Numerically calculated magnetostatic mode profiles in a cubic sample of edge size 40a; 
the profiles are depicted along the central axis (indicated in Fig. EJ) showing separately: (a) the 
dynamical magnetization m + relative values, and (b) the corresponding absolute values, \m + \. 



FIG. 5: The modes shown in Fig. are divided here into symmetrical (right column) and an- 
tisymmetrical (left column) groups. The modes within each group can be regarded as a family 
generated by the bulk- extended top mode: with lowering mode frequency, the ferromagnetic mode 
tails shrink, resulting in the appearance of three mode localization types: central-bulk, non-central 
and surface localization. 



FIG. 3: The quantized magnetostatic mode frequencies vs. the cube size, N; m indicates the mode 
number. 
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FIG. 6: The quantized magnetostatic mode frequencies (square points) vs. the mode number, m; 
note the double degeneration of the lower-spectrum modes. Note also that the frequencies of modes 
m = N and m = N — 1 are clearly separated from those of the remaining modes. The bold line 
represents the spatial distribution of the local dipolar field Q n values along the central axis. The 
cube edge size is (a) 2L = N - 1 =40, and (b) 2L = N - 1 =100. 



FIG. 7: Juxtaposition of magnetostatic mode profiles in two cubes of different sizes. Note that 
when the cube size increases, the lower modes tend to localize in the surface region, while the upper 
modes become strongly localized around the cube center; the two highest modes remain practically 
unchanged. 



FIG. 8: An illustration of the size-induced localization effect. When the cube size increases, the 
localization regions of the localized modes (m = N-Q, 5 and 1) narrow down substantially. 



FIG. 9: Magnetostatic mode spectrum in a cube of size 40a confronted with the spatial distribution 
of the local dipolar field O ra (bold line) along the cube central axis. A striking feature is that the 
local field line borders the "empty" region where mode amplitudes are suppressed. The bulk- 
extended mode frequencies lie above the local field range. 
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